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Abstract 

We explore an alternative strategy to determine the neutrino mass hierarchy by making use 
of possible future neutrino facilities at Fermilab. Here, we use CPT-conjugate neutrino channels, 
exploiting a Un beam from the NuMI beamline and a v e beam from a betabeam experimental setup. 
Both experiments are performed at approximately the same (E)/L. We present different possible 
accelerator scenarios for the betabeam neutrino setup and fluxes. This CPT-conjugate neutrino 
channel scenario can extract the neutrino mass hierarchy down to sin 2 29\3 ~ 0.02. 
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I. INTRODUCTION 



During the last several years thephysics of neutrinos has achieved remarkable progress. 
The experiments with solar , atmospheric {3], reactor and also long- 

baseline accelerator , [ll| neutrinos, have provided compelling evidence for the existence 
of neutrino oscillations, implying non zero neutrino masses. The present data require two 
large (#12 and #23) and one small (#13) angles in the neutrino mixing matrix Jjj], and at least 
two mass squared differences, Am] ; = m 2 — mf (where m/s are the neutrino masses), one 
driving the atmospheric (Am^) and the other one the solar (Am^) neutrino oscillations. 
The mixing angles 9 12 and #23 control the solar and the dominant atmospheric neutrino 
oscillations, while ^3 is the angle limited by the data from the CHOOZ and Palo Verde 



reactor experiments 

The Super-Kamiokande (SK) and K2K 9] data are well described in terms of dominant 
Vn — > v T ipp — * v T ) vacuum oscillations. A recent global fit [15| provides the following 3a 
allowed ranges for the atmospheric mixing parameters 

I Am§! I = (2 - 3.2) x 10" 3 eV 2 , 0.32 < sin 2 9 23 < 0.64 . (1) 

The sign of Am| 1; sign(Am| 1 ), cannot be determined with the existing data. The two 
possibilities, Amfj > or Am|j < 0, correspond to two different types of neutrino mass 
ordering: normal hierarchy and inverted hierarchy. In addition, information on the octant 
in which #23 lies, if sin 2 2#23 7^ 1, is beyond the reach of present experiments. 

The 2-neutrino oscillation analysis of the solar neutrino data, including the results from 
the complete salt phase of the Sudbury Neutrino Observatory (SNO) experiment 6j, in 
combination with the KamLAND spectrum data 1161, shows that the solar neutrino oscilla- 

n 

tion parameters lie in the low-LMA (Large Mixing Angle) region, with best fit values |15j 
Am 2 x = 7.9 x 10~ 5 eV 2 and sin 2 9 12 = 0.30. 

A combined 3-neutrino oscillation analysis of the solar, atmospheric, reactor and long- 
baseline neutrino data 15| constrains the third mixing angle to be sin 2 9 13 < 0.04 at the 3<r 
CL. However, the bound on sin 2 #13 is dependent on the precise value of Am^. 

The future goals for the study of neutrino properties is to precisely determine the already 
measured oscillation parameters and to obtain information on the unknown ones: namely 
#13, the CP-violating phase 5 and the type of neutrino mass hierarchy (or equivalently 
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sign(Am| 1 )). In the presence of matter effects, the neutrino (antineutrino) oscillation prob- 
ability gets enhanced [l?], for the normal (inverted) hierarchy. Making use of the different 
matter effects for neutrinos and antineutrinos seems, in principle, the most promising way 
to distinguish among the two possibilities: normal versus inverted hierarchy. However, the 
sensitivity to the mass hierarchy determination from the neutrino-antineutrino comparison 
is highly dependent on the value of the CP violating phase. Thus, possible alternative 
methods were first proposed in Ref. [l9]. In this paper we concentrate on the extraction of 
the neutrino mass hierarchy by combining a — > v e experiment with its CPT conjugated 
channel i> e — > i>^, see Ref. [19[. More recently, it is primarily the CPT-conjugate channel 



pairs that give the CERN-MEMPHYS proposal sensitivity to the hierarchy, see Ref. [20|. If 
nature respects CPT symmetry, then, at the same E/L the only difference between the two 
flavor transitions can come from matter effects and that near the first oscillation maximum 

P{y^ — > u e ) > P{v e — > Vp) for Normal Hierarchy 
and P{Vfj, — » v e ) < Pipe for Inverted Hierarchy, 

i.e. for the normal hierarchy the neutrino channel is enhanced and the antineutrino CPT 
conjugate channel suppressed and vice versa for the inverted hierarchy. This is the effect 
that will be exploited in this paper to determine the neutrino mass hierarchy. 

We will show that the combination of the Phase I (neutrino-data only) of the long-baseline 



v e appearance experiment 



NOt/A |2l|, e xploiting the off-axis technique 1 with a possible 



future betabeam facility 23j, |24], |25|, |26|, |27| at Fermilab exploiting a v e neutrino beam 



from radiative ion decays could help enormously in measuring the neutrino mass hierarchy. 
For our analysis, unless otherwise stated, we will use a representative value of lAmg-J = 
Z.o x ICT 3 eV 2 and sin 2# 23 = 1. For the solar oscillation parameters Am^ and #12, we will 
use the best fit values quoted earlier in this section. The structure of the paper is as follows. 
In Section [Til we present the general physics strategy used to determine the neutrino mass 



1 A neutrino beam with narrow energy spectrum can be produced by placing the detector off-axis, i. e., 
at some angle with respect to the forward direction. The resulting neutrino spectrum is very narrow in 
energy (nearly monochromatic, AE/E ~ 15 — 25%) and peaked at lower energies with respect to the 
on-axis one. The off-axis technique allows a discrimination between the peaked v e oscillation signal and 



the intrinsic v e background which has a broad energy spectrum [22]. In addition, the off-axis technique 
reduces significantly the background resulting from neutral current interactions of higher energy neutrinos 
with a 7T° in the final state. 
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hierarchy including the CPT conjugate channels used in this paper. Section HTT1 contains a 
realistic description of possible future betabeam facilities at Fermilab. The different scenarios 
deal with different ions, baselines and luminosities, and the performance of the strategy 
followed here in each of these scenarios is illustrated in Section IIVI The sensitivity curves 
for the several scenarios will be presented in Section |V] and the final remarks are summarized 
in Section IVII In the Appendix [Aj we discuss the details associated with comparing CPT 
conjugate neutrino oscillation probabilities. 



II. COMBINING NEUTRINO CHANNELS 



The strategy we have introduced in the previous section and we explain in detail here is 
different from the usual one, which exploits the combination of the neutrino and antineutrino 
oscillation channels. Typically, the proposed long baseline neutrino oscillation experiments 
have a single far detector and plan to run with the beam in two different modes, muon 
neutrinos and muon antineutrinos. In principle, by measuring the probability of neutrino and 
antineutrino flavor conversion, the values of the CP-violating phase 8 and the sign(Am| 1 ) 
could be extracted, since, in the presence of matter effects there will be two allowed regions 
for each type of hierarchy, normal or inverted, in the Piv^ — > v e ) versus Pip^ — > v e ) plane. 
In practice, the neutrino-antineutrino comparison does not provide the ideal tool to extract 
the neutrino mas hierarchy, as we explain below. 

Suppose we compute the oscillation probabilities Piy^ — > u e ) and P{p^ — > u e ) for a given 
set of oscillation parameters and the CP-phase 5 is varied between and 2ir: we obtain a 
closed CP trajectory (an ellipse) in the bi-probability space of neutrino and antineutrino 



conversion 



281 ] . Matter effects are responsible for the departure of the center of the ellipses 



from the diagonal line in the bi-probability plane for normal and inverted hierarchy. In 
Figure [D, we have illustrated the case for E = 2.0 GeV and L = 810 km, which roughly 
correspond to those of the NOi^A experiment [2lJ. The distance between the center of the 
ellipse for the normal hierarchy (lower blue) and that for the inverted hierarchy (upper red) 
is governed by the size of the matter effects. Notice that the ellipses overlap for a significant 
fraction of values of the CP-phase 5 for every allowed value of sin 2 29\3. This makes the 
determination of sig^Am^) extremely difficult, i. e., the sign (Am^) -extraction is not free 
of degeneracies. 
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FIG. 1: The bi-probability plot for P{v^ — > v e ) versus P(vu — ► & e ) at a baseline of 810 km and 
an energy of 2.0 GeV for the normal (blue) and the inverted (red) hierarchies. The smaller, lower 
(larger, upper) ellipses are for sin 2 2#i3 = 0.02 ( 0.10J. 



Following the line of thought developed by Minakata, Nunokawa and Parke [191] . we 
exploited in a previous work 29J, [30] the neutrino data only from two experiments at different 
distances and at different off-axis locations, such that the (E)/L is the same for the two 
experiments (see also Refs. [3l|, |32|, |33|, |3J, [35|)- In the case of bi-probability plots for 



neutrino-neutrino modes at different distances (which will be referred as near (N) and far 
(F)), the CP-trajectory is also elliptical. In Figure [2] (a) we present the bi-probability 
plot for the mean energies and baselines of the v e appearance experiments T2K {ssj and 



NOiM. [21[. The overlap of the two ellipses, which implies the presence of a degeneracy of 
the type of hierarchy with other parameters, is determined by their width and the difference 
in the slopes. Using the fact that matter effects are small (aL A 31 , being a = GpN e / a/2 w 
(4000 km)~ l the matter parameter), we can perform a perturbative expansion and assuming 
that the (E)/L of the near and far experiments is the same 2 , at first order, the ratio of the 
slopes reads [lJl 



aq 



1 + 4 (on-^n — apL F ) 



A. 



.31 



tan(A 



31, 



(2) 



The reason for this choice of (E) /L is explained in the next paragraph. 
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FIG. 2: (aj The left panel is the bi-probability plot for P(vu — > v e ) versus P{yu — > v e ) with baselines 
295 km and 810 km for the normal (blue) and the inverted (red) hierarchies. The smaller, lower 
(larger, upper) ellipses are for sin 2 28 '13 = 0.02 ( 0.10). The mean neutrino energies are chosen 
such that the {E) /L for the two experiments are approximately identical. 

(b ) The right panel is the bi-probability plot for P(vu — > v e ) versus P{p e — > Vy) for the normal (blue ) 
and the inverted (red) hierarchies. The baseline and mean neutrino energy for both experiments 
are 810 km and ~ 2 GeV, respectively. The smaller, lower (larger, upper) squashed ellipses are for 
sin 2 26»i3 = 0.02 (0.10). 

where a + and a_ are the slopes of the center of the ellipses as one varies #13 for normal 
and inverted hierarchies, ap and are the matter parameters, and Lp and are the 
baselines for the two experiments. The separation between the center of the ellipses for the 
two hierarchies increases as the difference in the matter parameter times the path length, 
(aL), for the two experiments increases. Also, since (A -1 — cot A) is a monotonically 
increasing function of A, we conclude that the smaller the energy, the larger the ratio 
of slopes, assuming the same (E)/L. However the width of the ellipses is crucial: even when 
the separation between the central axes of the two regions is substantial, if the ellipses for 
the normal and inverted hierarchy overlap, the hierarchy cannot be resolved for values of 
the CP phase, S, for which there is overlap. The width of the ellipses is determined by the 
difference in the (E) / L of the two experiments. 



6 



In the case of bi-probability plots for the — > v e and its CPT conjugated channel 
^ e — » &n at the same energy divided by baseline, (E)/L, the CPT-trajectory collapses to a 
line (see Figured (b)). As for the neutrino-neutrino case, we can perform a perturbative 
expansion, and, assuming that the (E)/L of the CPT conjugated channels is the same (to 
minimize the ellipses width), at first order, the ratio of the slopes reads (see Appendix and 
also Ref. Q) 

^±^l + A(aL + a CPT L CPT ) 1 ) , (3) 

a- V A 3i tan(A 31 )y 

where a + and a_ are the slopes of the center of the ellipses as one varies #13 for normal 
and inverted hierarchies, a and cicpt are the matter parameters and L and Lcpt are the 
baselines for the two experiments which exploit the — > v e and its CPT conjugated channel 
[p e — ■> v^). Notice that, compared to the neutrino-neutrino case given by Eq. (T5]), the 
separation between the center of the ellipses for the two hierarchies increases as the sum of 
the matter parameter times the baseline, aL, for both experiments does. Here the ratio of 
the slopes is enhanced by matter effects for both — > v e and its CPT conjugated channel 
v e — *• *V Figure [2](b) shows the bi-probability curves for the combination of these two flavor 
transitions, assuming that the two experiments are performed at the same mean energy and 
baseline. If the (E) /L of both experiments is the same, the ellipses will become lines with 
a negligible width. The separation of the lines for the normal and inverted hierarchy grows 
as the matter effects for both experiments increase. Consequently, the comparison of CPT 
conjugated channels is more sensitive to the neutrino mass hierarchy than the neutrino- 
neutrino one. 



III. BETA BEAMS AT FERMILAB 



A betabeam facility exploits a beam of electron neutrin os ( antineutrinos) from boosted- 
ion (3 + (/?") decays in the straight section of a storage ring [23j, |24j. The idea of considering 
higher 7 factors (and, consequently, longer detector baselines) was first proposed in Ref. 25]. 
An extensive phenomenological work has been devoted in order to optimize the betabeam 
physics reach, analyzing several scenarios with different 7 factors, boosted-ions and/or de- 
tector baselines [261 ] . 

Early on, 6 He and 18 Ne were identified as optimal ions, because of the low Q factor of 
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their decay. The lower the neutrino energy is in the rest frame, the more boost is needed 
to get to a given energy, and since the angular spread of the beam goes as I/7 this yields 
a more focused neutrino beam, which in turn produces more events in the far detector. 
Recently, it was proposed to use 8 Li and 8 B, which could potentially be produced in large 
amounts using a small storage ring with an internal gas target [381] . Since these ions have 
larger Q factor, they produce fewer neutrinos per ion in the far detector for a fixed neutrino 
energy and baseline. However, because less boost is needed, a smaller accelerator would 
be needed to achieve the same neutrino energy, as compared to the case of 6 He and 18 Ne. 
In this section we present possible betabeam scenarios at Fermilab, exploiting its current 
accelerator facilities. Since the analysis considered in this paper only requires anti-neutrinos 
from a beta-beam, we concentrate on 6 He and 8 Li. 

A rather thorough study of achievable ion intensities has been done at CERN [37]. The 
CERN scenario uses the existing PS and SPS accelerators, and would in addition require a 
proton source (e.g. the proposed Superconducting Proton Linac), target station, ion source, 
ion linac and Rapid Cycling Synchrotron (RCS), as well as a decay ring operating at SPS 
top energy. Based on a 6 He ion production rate of 2 x 10 13 /s, ions decaying in the straight 
section directed towards the experiment would produce approximately 3 x 10 18 antineutrinos 
per year. We consider two possible scenarios at Fermilab, namely: 1) accelerating 6 He to 
Tevatron top energy and 2) accelerating 8 Li to Main Injector (MI) top energy. These two 
scenarios produce neutrinos of comparable energies. The ions would be generated using a 
proton source (e.g. the Project X linac) and accelerated using e.g. a linac and a small RCS 
before being injected into the existing Booster. Possibly, the Recycler could also be used 
to accumulate bunches while the MI is ramping. In both cases, a new decay ring would be 
needed to store the ions 3 . 

Extrapolating from the work done at CERN, it appears reasonable to expect a useful ion 
decay rate (decays in the direction of the experiment) of about 1 x 10 18 6 He per year in the 
Tevatron case. At this intensity, the average power deposition in the Tevatron would be 
about 1 W/m, which is a generally accepted limit for hands-on-maintenance. Preliminary 
simulations indicate that the Tevatron magnets would be able to handle the distributed 
energy deposition from decay products in the arcs, but special care would have to be taken 



Note that if the Tevatron top energy is used ( 6 He), the decay ring would be very large and expensive. 
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to cope with the build-up of decay products in the straight sections. In the case of 8 Li in MI, 
the injectors could operate with a significantly higher duty factor, since there is no need to 
wait for the slow Tevatron ramp. However, at repetition rates and intensities corresponding 
to a useful ion decay rate of about 1 x 10 19 per year, activation of the Booster from decay 
products would become a serious issue. 

A very important property of the neutrino beam is the duty factor, defined as the relative 
fraction of time occupied by the neutrino pulse. This is used to suppress background by 
gating the data acquisition in the experiment. In the CERN study, a duty factor of a few 
per mil was obtained with considerable difficulty A small duty factor is challenging because 
it requires the ions to be concentrated in very few bunches, which among other things can 
cause space-charge problems, in particular at low energies. Using the CERN number of 
1 x 10 12 ions per RCS cycle, and a Booster injection energy of around 500 MeV/u, the beam 
must be distributed over about 10 bunches in the Booster to keep the space charge tune 
shift at an acceptable level. Approximately eight transfers from the booster per MI cycle 
would be required to obtain 1 x 10 19 useful 8 Li decays per year at MI top energy. Without 
RF manipulations, this would yield a neutrino beam duty factor of about 10%, but a duty 
factor of order 1% could likely be obtained by coalescing bunches at MI top energy 4 . In the 
case of 6 He, about twenty booster injections per cycle would be required per cycle, in order 
to compensate for the long Tevatron cycle and obtain a useful 6 He decay rate of 1 x 10 18 
per year. Assuming bunches are coalesced in the MI, this should also yield a neutrino beam 
duty factor of about 1%. Space charge is not expected to be an issue in the MI or Tevatron 
at these bunch intensities. 

For the 6 He case, therefore, it appears possible to generate 10 18 useful ion decays per year 
with a maximum gamma of 7n e = 350. For the 8 Li case, at a maximal gamma of 7u = 55, 
the rate could be higher (as explained above). We will explore an optimistic scenario of 
5 x 10 19 useful ion decays per year, as well a more pessimistic scenario of 10 19 useful ion 
decays per year from 8 Li decays. We will assume a duty factor of 1% for both ion species. 
Table [J shows the maximum Lorentz gamma factors in the Fermilab machines for the 6 He 
and 8 Li, as well as other ions considered in the literature. 

4 Although coalescing is standard procedure for generating single proton bunches from 53 MHz booster 
beam, the stability of the coalescing process needs to be demonstrated when generating multiple intense 
coalesced bunches simultaneously. 
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Machine 


6 He 2 + 


8 Li 3+ 


18 Ne 10+ 


8 B 5+ 


Linac 


1.6 


1.07 


1.15 


1.19 


Booster 


3.3 


3.7 


5.4 


6.1 


Main Injector 


54 


60 


90 


101 


Tevatron 


351 


395 


586 


659 



TABLE I: Maximum Lorentz gamma factors obtainable in the Fermilab machines for different ions 
of interest. 

IV. POSSIBLE EXPERIMENTAL SETUPS AT FERMILAB 

As we discussed in Section [Til the most sensitive, degeneracy free method, to extract 
the neutrino mass hierarchy exploiting a future high intensity conventional neutrino beam 
(i/p —>■ u e ) is the combination of this channel with its CPT conjugate u e — > v^. Future 
facilities like betabeams (neutrino factories) can produce neutrino beams which are entirely 
(partially) composed of v e or u e . 

The experimental strategy that we follow here is to combine the NOt/A experiment, which 
will measure the flavor transitions — > v e , with its CPT conjugated channel. The NOz/A 
experiment is expected to run at least five years with neutrinos. A 30 kton low density 
tracking calorimeter with an efficiency of 24% would be located at a baseline of 810 km and 
at 12 km off- axis distance from the beam center, resulting in a mean muon neutrino energy 
of 2.0 GeV. For the CPT conjugate channel, v e — > z/ M , we exploit possible, future betabeam 
facilities at Fermilab described in the previous section for two antineutrino emitters: 6 He and 
8 Li. A future neutrino factory exploiting neutrino fluxes from muon decays could also provide 
the v e CPT conjugate channel, if /i _ 's are stored in the decay ring. In the present study 
we do not explore this possibility. For the 6 He ion case, (maximum 7h c = 350), the mean 
electron antineutrino energy is ~ 1.2 GeV ((E u ) ~ 7He-^o,He, -Eo,He — 3.5 + m e MeV, being 
the electron end-point energy for 6 He). We present the results for two possible experimental 
6 He setups. In the first scenario, we consider a single baseline L = 810 km, at which the 
NOz/A detector would be located, with a total detector mass of 40 kton. This first scenario 
could be easily achieved adding to the 30 kton NOz/A far detector a second 10 kton detector 
at the NOz/A far site. Possible detector technologies are Liquid Argon or iron calorimeter 
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detectors. Liquid Argon (LAr) detectors have excellent efficiency, background rejection and 
energy resolution, but they could suffer from a large atmospheric neutrino background, which 
could be overcome only if the beam duty cycle is < 10 -4 24j. In the second scenario, we 
consider a detector similar to the one of the MINOS [10] experiment (5 kton) at 735 km but 
twice in size. If the ion luminosity could be improved by a factor of two, MINOS far detector 
would be sufficient. This scenario benefits from the lower atmospheric neutrino backgrounds 
at the MINOS site. The beam duty cycle would not, therefore, be a major issue, and a duty 
factor of ~ 1% would be sufficient to overcome the atmospheric neutrino background in this 
case. For the 8 Li ion case, (maximum '-fu = 50), the mean electron antineutrino energy is 
~ 1.8 GeV ((E u ) ~ 7Li-E'o,Li, E 0jLi = 16 MeV, being the electron end-point energy for 8 Li). In 
order to ensure an almost degeneracy free hierarchy measurement, the (E)/L of the — > v e 
channel from NOi^A and its CPT conjugate channel should be similar, therefore the detector 
should be located at L = 300 km. We will consider 5 x 10 19 (1 x 10 19 ) useful ion decays per 
year with a 10 (50) kton detector respectively. As previously discussed, considering 5 x 10 19 
useful ion decays is a quite optimistic assumption. However, notice that the same statistics 
could be achieved with the more conservative assumption of 1 x 10 19 useful ion decays per 
year, with a larger 50 kton detector (these two configurations provide the same statistics). 
Due to the lower neutrino energies, the detector could be a Liquid Argon TPC, a NOz/A-like 
Totally Active Scintillator Detector (TASD), or a water Cherenkov. For the NUMI off-axis 
neutrino beam, we have not considered a binning in the signal, since neutrino events will 
lie in a very narrow energy window [1.5,2.5] GeV. The 6 He and 8 Li antineutrino betabeams 
have been divided in three energy bins, assuming an energy detection threshold of 0.5 GeV. 
For the 6 He case, we divide the signal in three energy bins, in the energy ranges [0.5,1.0], 
[1.0,1.5] and [1.5,2.4] GeV, respectively. In the 8 Li case, the energy ranges are [0.5,0.75], 
[0.75,1.0] and [1.0,1.5] GeV, respectively. 

Figure [3] depict the bi-event plots for the combination of the NOz/A neutrino events 
(z/ M — > v e ) with its CPT conjugated channel (z/ e — > v^) from the betabeam experiment, 
resulting from the decays of 6 He (left panel) and 8 Li (right panel), for both normal and 
inverted hierarchies. The statistics considered for NOz/A correspond to Phase I of the 
experiment (neutrino running only). For the 6 He betabeam experiment, we assume a number 
of useful ion decays of 10 18 per year, five years of data taking, and a 40 kton detector located 
at 810 km. For the 8 Li betabeam experiment, we assume 5 x 10 19 (1 x 10 19 ) useful ion decays 
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FIG. 3: (a) The allowed regions in the bi-event plot for N{i>u — ► u e ) for NOvA versus N(v e — ► £v) 
/or a betabeam experiment which exploits antineutrinos from e He decays with 7# e = 350, and a 
detector of 40 kton located at a distance of L = 810 km. The blue (red) ellipses denote normal 
(inverted) hierarchies. From bottom up, the ellipses correspond to sin 2 29^ varying from 0.01 to 
0.1. The solid (dashed) ellipses illustrate the third (second) energy bins of the betabeam spectrum, 
(b) Same as (a) but with antineutrino fluxes resulting from 8 Li decays, and a detector of 10 kton 
at 300 km. 

per year, ten years of data taking, and a 10 (50) kton detector located at 300 km. 

Figure [31 (a) shows that, for the combination of NOiM. off-axis neutrino events with the 
antineutrino events from 6 He decays, the separation between the bi-event contours for the 
normal and inverted hierarchies is larger than in the case of 8 Li generated antineutrino 
events, as seen in Figure [3] (b). As previously explained, the difference in the slopes of the 
two hierarchies is proportional to the sum of the size of matter effects times the baseline, 
o-novaLnouA + o-cptLcpt- The product cicptLcpt is larger for the 6 He betabeam v e events 
(with a baseline of 810 km), than for the 8 Li betabeam v e events (with a baseline of 300 km). 
The solid (dashed) contours in Figure [3] show the number of betabeam antineutrino events 
in the second (third) energy bin. When the (E)/L of the —>■ v e and its CPT conjugated 
channel are similar, the ellipses width is minimal (they collapse to a line) and therefore the 
elliptical contours for the normal and inverted hierarchies will not overlap, regardless the 
value of the CP violating phase 5. For the combination of NOz/A off-axis neutrino events 
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with the 6 He betabeam antineutrino events, there exists a clear difference between the second 
and third energy bins in the bi-event contours: while they are ellipses in the former, they 
are almost lines in the latter. Only in the third energy bin ([1.5,2.5] GeV) is the (E)/L the 
same for the 6 He betabeam v e and for the NOz/A events. For the 8 Li case the ellipses 
width is minimal for both the second and third energy bins, since both bins are close to 
E ~ 0.8 GeV, the energy at which the ((E) jV)u equals the ((E) / L)novA- 

V. SIGN Amfi SENSITIVITIES 

In this section we present the physics results from the combination of antineutrino data, 
resulting from several possible betabeam setups, with the NOz/A neutrino data. 

Figure H] (a) shows the 90% C.L mass hierarchy sensitivity, assuming two degrees of 
freedom statistics (2 d.o.f, that is, A% 2 > 4.21) for the combination of NOz/A neutrino 
data with the 6 He betabeam antineutrino data, neglecting the background in the z/ e — > 

channel. In order to safely neglect the atmospheric neutrino background, one would 
need a very small, experimentally challenging duty factor for the betabeam neutrino fluxes. 
The baseline for the two experiments is fixed at 810 km and the 7n e = 350, in order to 
have a similar (E)/L in both the muon neutrino and the electron antineutrino channels. 
The binning of the signal has been chosen as quoted in the previous section. We have 
considered a flux of 10 18 antineutrinos per year, and five years of data taking. The blue 
(red) lines assume normal (inverted) hierarchy, and the solid (dotted) lines depict the results 
for a betabeam antineutrino experiment with a 40 (10) kton detector. For the 40 kton 
detector, the sensitivity is better for the inverted mass hierarchy, since in this case statistics 
is dominated by the antineutrino channel z/ e — > v^. For the 10 kton detector, the two 
channels (i.e., the neutrino channel from NOz/A and the antineutrino channel from the 
betabeam experiment) will have similar statistics and therefore the sensitivity for the normal 
and the inverted mass hierarchies are similar. As a comparison, we also show the results 
for the NOz/A experiment (upgraded by a factor of five in statistics), assuming five years 
of neutrino and five years of antineutrino data taking, see the dashed curves in Figure HI 
The setup proposed here improves the sensitivity of the NOz/A upgraded experiment by an 
order of magnitude, and more importantly, eliminates the dependence of the mass hierarchy 
determination on the value of the CP violating phase S. 
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Figure H] (b) shows the 90% CL mass hierarchy sensitivity, assuming two degrees of free- 
dom statistics (2 d.o.f, that is, A% 2 > 4.21) for the combination of NOz/A neutrino data 
with the 6 Li betabeam antineutrino data, neglecting the background in the v e — > z> M chan- 
nel. The baseline for the betabeam experiment is L = 300 km and the binning of the signal 
has been chosen as described in the previous section. We assume 5 x 10 19 8 Li generated 
antineutrinos per year, and ten years of data taking. The blue (red) lines assume normal 
(inverted) hierarchy, and the the solid (dotted) lines depict the results for a betabeam an- 
tineutrino experiment with a 10 (2) kton detector. For the case of 8 Li betabeam experiment, 
the sensitivity is similar for the normal and inverted mass hierarchy, and smaller than in the 
case of 6 He betabeam experiment. This is expected, since at a shorter baseline (300 km) 
the product of the matter potential times the distance is reduced. Again, the combination 
of the NOz/A neutrino data only with the 8 Li betabeam antineutrino data provides a much 
better sensitivity to the mass hierarchy than the NOz/A upgraded experiment alone. 

However, as previously stated, the beam duty cycle needed in order to neglect the at- 
mospheric neutrino background is highly challenging. For a MINOS-like detector, there are 
30 atmospheric neutrino interactions per kton-year which could mimic a muon coming from 
the oscillated v e —>■ [39]. In order to avoid such a large background, we have assumed 
a betabeam duty cycle ~ 10~ 2 , which seems experimentally achievable. Figure [5] shows 
the equivalent to Figure H] when adding the atmospheric neutrino background quoted above, 
rescaled accordingly to the detector sizes and the exposure times. Notice that the presence 
of a non negligible atmospheric neutrino background in the antineutrino channel v e — > 
reduces the sensitivity reach especially for the case of the inverted mass hierarchy. However, 
even in the presence of a non negligible background, and with the very conservative assump- 
tion of a 10~ 2 beam duty cycle, the combination of the betabeam antineutrino data with the 
NOz/A Phase I neutrino data, provides better sensitivity than the NOz/A experiment alone 
(upgraded by a factor of five and running in both neutrino and antineutrino mode). If a 
smaller duty factor < 10~ 2 could be achievable (as commonly assumed, following Ref. 24j ) 
the sensitivity to the mass hierarchy would lie within the limits illustrated in Figure H] (the 
most optimistic case with no atmospheric neutrino induced background) and the limits de- 
picted in Figure [5] (the most pessimistic case with atmpospheric neutrino backgrounds and 
a beam duty cycle ~ 10~ 2 ). 
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FIG. 4: (a) The 90% CL (2 d.o.f) hierarchy resolution curves for different exposures for the ®He 
betabeam v e fluxes at 810 km, combined with five years of neutrino data only from the NOvA far 
detector, located 12 km off-axis at 810 km. Only backgrounds in the NOvA experiment have been 
included. The blue (red) curves assume normal (inverted) hierarchy. The solid (dotted) line depicts 
the results for 2 x 10 20 (5 x W 19 ) useful ion decays times kton. The blue (red) dashed curve shows 
the sensitivity reach at 90% CL (2 d.o.f.) from the combination of neutrino and antineutrino data 
from the NOvA experiment, assuming five years running in both neutrinos and antineutrinos with 
a factor of five increase in statistics and normal (inverted) hierarchy, (b) The 90% CL (2 d.o.f) 
hierarchy resolution curves for different exposures for the 8 Li betabeam electron antineutrino fluxes 
at 300 km, combined with muon neutrino data only from the NOvA far detector at 12 km off-axis 
at 810 km. Again, only backgrounds in the NOvA experiment have been included. The blue (red) 
curves assume normal (inverted) hierarchy. The solid (dotted) line depicts the results for 5 x 10 21 
(10 21 ) useful ion decays times kton. 

VI. CONCLUSIONS 

We have explored an alternative strategy for measuring the neutrino mass hierarchy. 
Unlike the approach followed by future long baseline neutrino oscillation experiments that 
combine the neutrino-antineutrino data, the combination of the CPT conjugated channels 
that we study here provides an almost degeneracy free determination of the neutrino mass hi- 
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FIG. 5: Same as Figs. [7] but including backgrounds in the betabeam electron antineutrino data, see 
text for details. 

erarchy, provided the two channels have similar (E)/L. Future neutrino facilities at Fermilab 
could provide these CPT neutrino-conjugated channels. The NOiM. v e off-axis appearance 
experiment could provide the —>■ u e channel. A future betabeam facility based at Fermilab 
could provide the CPT-conjugated v e — > channel. A realistic estimate of the expected 
electron antineutrino fluxes from boosted ion decays is presented. We propose two possible 
accelerator scenarios for generating the betabeam electron antineutrino fluxes: the Tevatron, 
which could accelerate 6 He ions, and the Main Injector, which could accelerate 8 Li ions. In 
the case of the Tevatron, the decay ring would be very large and possibly prohibitively 
expensive. The first scenario could benefit from the NOz/A far detector at L = 810 km, 
(but the decay ring needed would be very large) ; for the second scenario, an additional, 
although smaller 2 — 10 kton MINOS like detector at a shorter baseline, L = 300 km, would 
be necessary (the decay ring needed in this case would be smaller, though). In the more 
pessimistic case, with a modest beam duty cycle of 10~ 2 and including realistic atmospheric 
neutrino backgrounds, the neutrino mass hierarchy could be determined for sin 2 2# 13 > 0.01, 
independently of the value of the CP violating phase 5, for both accelerator possibilities. 
These two alternative choices could improved by an order of magnitude the sensitivity to 
the neutrino mass hierarchy obtained by a future NOz/A upgraded experiment exploiting 
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both neutrinos and antineutrinos. 
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APPENDIX A: CPT CONJUGATE PROBABILITIES: 

The amplitudes for — > u e and u e — > consists of two terms, one associated with 
the atmospheric 5m 2 scale and the other associated with the solar 5m 2 scale. Thus, the 
probability for these the CPT-conjugate processes contain three terms; the square of each of 
the amplitudes plus the interference term between the two amplitudes which depends on the 
CP phase 5. For the normal (upper sign) and inverted (lower sign) hierarchy, the — > v e 
and u e — > appearance probabilities are given by 



where Ay = \5m 2 j\L/AE and a = GfN £ /\/2 (4000 km) -1 . The atmospheric amplitude 
for — > v e is ±^X±9 whereas the solar amplitude is \^Pq and the relative phase between 
these two amplitudes 5 is (±A 31 + 5). In vacuum, X + = X_ = X and the two probabilities 
are identical, as they must since they are CPT conjugates. 

5 The full amplitude for ^ -> v e is (±^/X^9 e -i(±A 3 i+«) + ^/p~). 




(Al) 



The coefficients P Q and X± are simply 
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The other related CPT conjugate pair of appearance probabilities, P[y e — > v^) and 
Pip 11 v e ), can be obtained from the above by changing the sign of S, as follows 

P(v e -> ^) = X ± 6 2 ± 2v^Iv^ ^cos(±A 31 - 5) + P G 





P(^ e -> ^) = X T # 2 ± 2v^v^ ^cos(±A 31 - 5) + P . (A2) 
The difference between the first two CPT conjugate appearance probabilities, is given by 



P{V» - ^e) " Pip e - ^) = ±9 (y/X^ 



X + + VX_)^ ± 2^/P Q cos(±A 13 + 5) 



This quantity is positive for the normal hierarchy (NH) and negative for the inverted hier- 
archy (IH), if 

(A3) 



(A4) 



VX+ > y/X. and 9 > 2^P Q /(y/X + + V x -) ~ V^o/V^o, 

for all values of the CP phase 5. The constraint on 9 requires 6 

sin 2 29 12 A 2 21 



sin 2 2# 13 > 



0.001 - 0.002, 



tan 2 #23 sin 2 A 3 i 

whereas the constraint, y/X + > y/X_, is satisfied near the first oscillation maximum pro- 
vided (aL) < 1, i.e. L < 4000 km. 

With these rather weak constraints then 



Piv„ 



u„) > P(u e 



for NH 



and P(i/ M -> z/ e ) < P(P e -> P M ) for IH 



(A5) 
(A6) 



for all values of the CP phase 8. For the normal (inverted) hierarchy, the matter effect 
enhances (suppresses) the Piy^ — > v e ) channel and suppresses (enhances) the Pip e — > v^) 
channel, thus the matter effect in a sense is used twice. Of course, the difference between 
these two appearance probabilities is larger at larger values of 9 and at larger values of 
the matter effect. This is the effect that is exploited here to determine the neutrino mass 
hierarchy. 

In the Piy^ — > u e ) versus Pip e — > D^) plane the trajectory for fixed value of 9 as the CP 
phase 5 is varied from to 27r is in general an ellipse which collapses to a line if the E/L of 
both channels is the same. The centre of this ellipse is given by 



(P(z/ e - Vj, P(^ - v e )) = (X T 9 2 + P , X ± 9 2 + P ). 



(A7) 



This is the value of 8 at which the atmospheric and solar amplitudes have the same magnitude in vacuum. 
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Thus, as 9 is varied, the centre of the ellipses form lines with slope given by 

a + = — forNH 
and a_ = ^± for IH. (A8) 

VS. _ 

If the matter effect is small, (aL) <C A 31 , one can perform a Taylor series about the vacuum 
such that 

ct ± = 1 =F 4(aL) [A~l - cot A 31 ] + O(aL) 2 . (A9) 

It is the difference in the slopes of the two lines (for the normal hierarchy, a + and for the 
inverted hierarchy, a_) which provides the separation between the allowed regions for two 
hierarchies in the Piy^ — > u e ) versus P{p e — > v^) plane. 
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